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We develop an interpretation of the off-equilibrium dynamical solution of mean-field glassy mod- 
els in terms of quasi-equilibrium concepts. We show that the relaxation of the "thermoremanent 
magnetization" follows a generalized version of the Onsager regression postulate of induced fluc- 
tuations. We then find the rationale for the equality between the fluctuation-dissipation ratio and 
the rate of growth of the configurational entropy close to the asymptotic state, found empirically in 
■ mean-field solutions. 
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Ch ' I. INTRODUCTION 

c ' 

Aging is a scaling dynamical regime characteristic of glassy systems [Q. In this regime, typical features of equilibrium 
systems, such as the asymptotic absence of macroscopic heat currents, coexist with non-stationary aspects such as 
, t | the dependence of the correlation and response functions on the system's "age", i.e. the time spent in the low 
temperature phase. The solution of mean-field spin glass models has given a general framework to understand 
aging phenomena, and has produced detailed predictions, which have been verified in numerical simulations of long- 
range |(| and short-range [?J glassy systems. A characteristic prediction of this solution is the existence, at low 
temperature, of a dynamical regime where extensive quantities depending on the configuration of the system at a 
single time (one time observables in the following) are well thermalized, (or evolve very slowly), while two time 
£j ■ correlation function and susceptibilities exhibit non-stationary scale invariant behavior. 
— 1 \ Despite this coherent theoretical scheme, and recent progress in linking aging dynamics to the nature of the equi- 
librium regime Q, a physical understanding of some fundamental aspects of aging dynamics is still lacking. While 
the investigation of both equilibrium and asymptotic off-equilibrium regimes give solutions that show unexpected 
qq coincidences, all efforts to interpret aging as a quasi-equilibrium condition have been thwarted by facts such as: 

' 1. No matter how large we take the aging time if we then wait long enough the system eventually wanders away 

, from any finite region of phase space |2j . 

O ■ 

2. Two identical systems starting from the same condition at any given aging time will always come apart as far 
Q\ \ as possible |J . 

In this paper we try to make ends meet. In particular, we try to give a physical and intuitive explanation of the link 
between the so called "fluctuation dissipation ratio" (FDR) (the factor x(q)) and the Parisi function (in Sherrington- 
Kirkpartick like models) or the derivative of the configurational entropy close to the threshold state (in p-spin-like 
models). Our work should be understood as a physical interpretation of mean field aging dynamics in terms of a 
Ch , quasi-equilibrium scenario and not as an alternative derivation of the results of the theory. 

We show that a modified version of the Onsager postulate on regression of fluctuations applies to the relaxation of 
the magnetization in thermoremanent magnetization experiments. The observed anomalies in the response are then 
analyzed. In this paper we discuss aging mainly in the case of a "one time sector" approximation, which is the dynamic 
counterpart of the "one step replica symmetry breaking" (1RSB) approximation in the equilibrium analysis. This is 
exact in models like the p-spin model, while it is only an approximation, and a rather crude one, when continuous 
replica symmetry breaking is present, like in the Sherrington-Kirkpatrick (SK) model or the random manifold model. 
We will refer to the first class of models as "p-spin-like" and to the second as "SK-like" . We discuss only this "one 
time sector" approximation to unify the argument, and simplify the notation. But we expect the reader to be able to 
generalize it without effort. 

The picture that emerges from our analysis is simple and intuitive: the age of an aging systems determines the 
rate of entropy decrease, i.e. the flow rate of heat towards the thermal bath. A small force in the linear response 
regime cannot change this rate. Acquiring a non zero magnetization means entropy reduction which has then to 
be compensated by an increase (or reduced decrease) of the free-energy associated with the spin-couplings. As a 
consequence, the response becomes proportional to the growth of the logarithm of the number of quasi-states (to be 
defined later) with free-energy. It is as if the slow degrees of freedom respond to external forces by sampling states 
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that lie above in free-energy while they are blocked from exploring those that are at the same or at a lower level. 
This paper presents what we believe are convincing arguments in favor of this assertion. As a byproduct the time 
scale dependent effective temperatures will appear jlO| and their connection with the derivative of the logarithm of 
the number of quasi-states (or configurational entropy) with respect to free-energy is explained. |Tl| l 

We organize the paper as follows: in section II we review some properties of aging in mean-field. In section III we 
discuss our definition of quasi-equilibrium and how it relates to dynamics. In section IV we recall Onsager regression 
postulate and generalize it to aging systems. In section V we discuss the origin of the FDR in SK-like models, while 
we treat the case of p-spin-like models in section VI. Finally we present a summary and the conclusions. 



II. A SHORT REVIEW 



We consider a mean-field spin-glass quenched at a given time t = into the glassy phase. For simplicity we will 
imagine that the degrees of freedom consist of Ising spins Si — ±1, i — 1,...,N (N — > oo). We are interested on 
the long time dynamics of such a system, i.e. on a regime where it has thermalized for a long time t w before any 
measurements. The long time limit t w — > oo is always taken after the thermodynamic limit. 

The free-energy, and its derivatives with respect to the control parameters (e.g. temperature) tend to some asymp- 
totic time-independent values. More interesting is the behavior of observables depending on two time variables, such 
as correlation and response functions. Aging behavior manifests itself as an asymptotic non-stationary dependence 
on time of these quantities. 

Let us consider the spin-spin time dependent autocorrelation function 

1 N 

C(t,u) = — Y,Si(t)Si(u) t>u»t w (1) 

i=l 

A first, short time, dynamical regime is obtained considering the difference r = t — u finite to derive stationary 
correlations C* t (r). Let us denote qEA the long r limit of C* t . In the aging regime C(t, u) relaxes below qEA- In the 
"one time sector" approximation scheme, the correlation function in this regime can be written as , 

C(t,u) = C ag (h(u)/h(t)) (2) 

where h(-) is an increasing function not derivable from the present theory, and t, u are large with lim^u—joo h(u) /h(t) = 
finite. The formulae in the two regimes are summarized in: 

C(t, u) = C st (t -u) + C ag (h(u)/h(t)) (3) 

where C s t{t — u) — C* t (t — u) — qEA is a monotonically decreasing function equal to 1 — qEA for t — u = and tending to 
zero for t — u — > oo, while C ag (h(u) /h(<)) is equal to qEA for h(u)/h(t) = 1 (u = t) and tends to go for n (u)/h(£) — > 0. 
In order to simplify the notation we will suppose that the "time reparametrization" h(i) is the identity h(t) = t, and 
C ag (h(u)/b.(t)) = C ag {u/t). We will also suppose go = 0, but this will not affect any of our following arguments. 

We stress that the form (||) of the correlation function implies that if we fix the value of u and let t run, the time 
spent at the value of the correlation equal to qEA is much larger then the time needed to reach it. This is an aspect 
of the "time scale separation" observed in glassy systems that will play a crucial role in our discussion. 

We are also interested in the behavior of the linear response function, 

I N 5 

^'^iEAm^^l" t>u»t w (4) 
N f-f 5hi(u) 

and the corresponding integrated function 

x(t,u) = J^ dsR(t,s) t>u»t w (5) 

which represent the susceptibility at time t in a "thermoremanent magnetization" experiment in which a constant 
small field has acted from the time t w up to time u. By the condition u » t w we mean t w /u — > for t w — > oo. The 
linear response theory requires that the limit h — * be taken before sending the time u to infinity. 

While in the stationary regime the fluctuation dissipation relation R(t, s) = R s t(t — s) — /3 8C g* :S ^ is verified, in the 
aging regime the relation is substituted by having a non trivial fluctuation-dissipation ratio (FDR): 
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x ^ = ,,1™. ^§MT ' (6) 

C(t,s) = q $ s 

which turns out to coincide with the function x(q) that appears in the replica approach that in principle applies to 
equilibrium of different kinds Jl2| , p^[ . In the one sector scenario x(q) is constant all through the aging regime. 

The FDR x(q) verifies the mathematical properties of a cumulative probability distribution, a feature that has been 
explained in recent work where it has been shown that there is a deep connection between the dynamic properties 
during aging and the property of ergodicity breaking in equilibrium |g] . Using only the hypothesis of equilibration 
of one time observables (OTO) and the existence of a linear response regime for the correlation functions (stochastic 
stability), it was proved that the function x(q) is related to the function P(q) describing the statistics of equilibrium 
pure states Q through the equation: 

Pi,) = (7) 
aq 

The theorem was originally formulated for finite dimensional systems, where the OTOs are guaranteed to thermalize 
but can be generalized to mean-field long range models of the SK-like class. A different situation is found in p-spin- 



like models |15 where one of the hypotheses of the theorem is violated because the asymptotic value of the dynamic 



energy is higher than the one of the states dominating the partition function. 



III. THE QUASI-EQUILIBRIUM HYPOTHESIS 

From now on we will work on the "one time sector" approximation described in the previous section. 

Let us consider a large time u, and the corresponding spin configuration Si(u). Our previous observations suggest 
that the value of qEA can be used to decompose the spin configuration in a "fast part" and a "slow part" according 
to 

Si(u) = [Si(u) - rrii(u)] + m^u) (8) 
where the slow variable rrii(u) can be estimated immediately from the running average 

j J-U+At 

m,i(u) ~ — / dw Si(w) with C(u + At, u) = qEA- (9) 

^ Ju 

We will see later how to improve on this estimate. By (||) and ([)]) we obtain correctly the corresponding decomposition 
of the correlation function in two time domains 

([Si(t) - mi(t)][S((u) - rmiu)}) = C st (t - u) (10) 

while 

(mi(t)mi(u)) = C ag (u/t). (11) 

any other decomposition, obtained averaging the spins over times such that the value of C(it, v) is different from qEA, 
would mix C s t and C ag - 

Through this decomposition one can define a dynamical notion of "quasi-state" in which the system (almost) 
equilibrates before relaxing further. The quasi equilibrium hypothesis can be formulated considering the probability 
distribution of finding the system in a given configuration of the slow and the fast variables at time t induced by the 
thermal noise and the flat distribution of initial conditions 

P t ({S i },{m i }) = {T\6(S i (t)-S i )6(m i (t)-rn i )) t ,™« . (12) 

A initial conditions 

i 

which can be written as: 

Pt({S i },{m i }) = P t ({S i }\{m i }))Pt({m i })). (13) 
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All the known properties of the dynamical solution, and in particular the short time response to external perturbations, 
are consistent with the proposition that the conditional probability P t ({Si}\{mi}) becomes independent of time, and 
takes asymptotically the form of a restricted Boltzmann measure^ 

We stress that the same property manifestly would not hold had we chosen time scales such that C(t, u) < Qea in 
the average (^). This proposition is implicit in the definition of the aging regime, where there is no ambiguity in 
the definition of qEA- The measure should be restricted to the transverse configuration space projecting out those 
directions along which the system evolves (transverse quasi-states) where the energy landscape is flat or has negative 
eigenvalues. In these conditions, the free-energy of the quasistate (see below), as well as the value of qEA entering in 
( pT| ) are close to their asymptotic values, but still depend on time. We can safely assume that in the asymptotic regime 
the number of negative directions become vanishingly small. Notice that, if we take two macroscopically different 
sets of slow variables {rrii} and {m^} then, by construction, the corresponding conditional probabilities P ({<Si}|{ m i}) 
and P ({<Si}|{m^}) are mutually orthogonal. This can be easily understood from the fact that the mutual overlap 
among a configuration with non-zero weight in the first distribution and a configuration with non-zero weight in the 
second one is almost surely smaller then qEA- It is therefore convenient to think of a discretized mi-sphere such that 
the centers of the neighboring cells correspond to disjoint quasistates. In general, we expect different quasistates to 
define disjoint regions in configuration space and that the union of all such regions define a partition of all relevant 
configurations. We will use a as the index of the quasistate which of course will change with the slow time. In (||) 
therefore rrii(u) should rather read mf with a a function of u. The finiteness of At limits the accuracy of the running 
average estimate. To derive a better estimate of mf we could clone the trajectory from time u on and take a weighted 
average along all trajectories. 

It is useful to define thermodynamic quantities such as the dynamical free-energy: 

?t = Y, P *({<D [F({mf})+T log (P t ({mf}))} (15) 

a 

where T is the temperature of the thermal bath, and 

Him?}) =/lI^ P ({^}K<» W{Si}) + TlogiP({Si}\{<}))] (16) 

i 

is the free-energy of the (transverse) a quasi-state. 

Observe that Tt includes the average free-energy of the quasi-states 

F t = £ P t ({mf})F({mf}) (17) 

a 

and a slow entropy term 

St = -Y. p t({<}) (Pt({mf})) . (18) 



An explicit computation shows that, due to the disjointness property of the quasistates, the sum St + St coincides 
with the entropy of the distribution P t ({Si}). 

We can identify T t with the dynamical free-energy / fli dSiP t ({Si})[H({Si}) + T log P t ({$})]. 

This last quantity is known to decrease in any process verifying detailed balance. In our case, due to the white 
average over the initial conditions we expect in addition both F t and St to decrease with time. 

For typical trajectories extensive quantities are self-averaging and therefore the free-energy F t is a well defined 
function. The asymptotic value, i 7 ^, is the free-energy of the equilibrium state for SK-like systems, and the free- 
energy of the threshold TAP solutions for the p-spin class. 



x To be more precise we should define the measure in such a way that each Si has average m;. One can check with a detailed 
calculation that this condition is automatically fulfilled by the measure (p^). 
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The role of St in the dynamical relaxation is not immediate. By construction, it does not say anything about the 
number of quasi-states accessible starting from a generic point of a trajectory at time t. In fact we expect that the 
support of Pt({rrii}) decomposes in non-overlapping, mutually inaccessible, regions of phase space that become more 
and more isolated as time advances. 

By inverting the relation F t among free-energy and time, we can define S(F) — St(F) an d derive that S(F oc ) 
is the extensive part of the configurational entropy of ground states in SK-like models (zero in this case) and the 
configurational entropy of the threshold states in the p-spin. The dynamical entropy that we defined weights different 
regions of phase space according to their basins of attraction. We are here using the observation that being the 
threshold states identical in all their properties they must also have the same basin of attraction. Furthermore with 
respect to those states which appear in exponentially smaller numbers there is the additional observation that their 
basin of attraction cannot be so much larger as to compensate for their smaller multiplicity. 

For large time, the support of Pt({rrii}) will be in regions of small TAP gradient. We can calculate S(F): 

St = S(F t ) = ^p-\ Foo (Ft - Foe) + Sco (19) 

In our scenario S (Ft) measures the multiplicity of quasistates at t which as before will have equal basin of attraction. 
A detailed calculation in the p-spin model is developed in the Appendix. It shows that if we compute the multiplicity 
of minima of the modulus of the gradient of the TAP free energy then their number and derivative are continuous 
across the threshold value. Therefore we propose to identify the quasistates with the minima of the gradient of the 
TAP free energy. For SK-like models we conjecture instead S(F) to be equal to the number of stable TAP excites 
states at level F. 

With this identification we now write 9 gp \ equal to (3x. It is one of those remarkable coincidences, referred to 
in the introduction, that the same value of x appears in the anomalous FDR. We will interpret this coincidence in 
the following sections. 

This identification of the quasistates is crucial. An explicit calculation of the dynamical entropy could check its 
validity but unfortunately with our present techniques such calculation is not feasible. 

A strong hint in favor of it comes however from the study of the equilibrium dynamics at temperatures slightly 
larger then the dynamical transition temperature Td- For T — Td « Td there is a similar separation of 2 time scales 
controlled by the ratio T T r Td . This allows the definition of dynamical quasistates along lines similar to the ones followed 
in the non-equilibrium case. The free-energy obtained by considering the collection of the quasistates of appropriate 
energy, should be coherent with the direct computation from the partition function. By an explicit computation that 
we sketch in the appendix, we verified that, up to second order corrections in T — Td, the thermodynamic entropy 
coincides with the TAP internal entropy plus the configurational entropy of the saddles. 

IV. REGRESSION OF FLUCTUATIONS AND ONSAGER POSTULATE 

In order to discuss the behavior of the response function we consider the set-up of "thermoremanent magnetization" 
(TRM) experiment jlj]. The system is allowed to age in a small field h acting from time t w to time u such that 
C(u,t w ) — > 0. At later times t > u one detects the magnetization M(t) = jfJ2iLi^i(^) = hx(t,u). ^ Our set-up 
differs slightly from the one considered usually in the literature, in which the field acts directly from the quenching time. 
We switch the field on at time t w because we find conceptually clearer to discuss the behavior of the magnetization 
starting from a situation where the system is already in the asymptotic regime. We notice that the response to any 
arbitrarily varying field hit) can be expressed as a linear superposition of TRM magnetizations. 

In order to discuss the decay of M(t) we will show that a generalization of the Onsager postulate of normal regression 
of fluctuations applies to the dynamical off-equilibrium process P6| . The principle, originally stated for equilibrium 
systems, concerns the behavior of macroscopic quantities and states that in the linear response regime one can not 
distinguish the regression of a spontaneous fluctuation of a certain quantity from the regression from the same value 
when imposed through a constraint on the equilibrium measure. Onsager's postulate means that for a large system, 
a spontaneous fluctuation must have the characteristic of the most probable fluctuations and therefore correspond to 
constrained minimization of the free-energy. This is equivalent to free-energy minimization in a conjugated field, thus 
leading to an immediate derivation of the fluctuation-dissipation theorem. 



2 It should be kept in mind that M(t) denotes the magnetization at time t but can depend on both t and u. 
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More in detail the argument can be phrased as follows . Consider a thermodynamic system at equilibrium and a 
given macroscopic (extensive) quantity a which takes the value zero at equilibrium. Be 7 the corresponding conjugate 
intensive variable. Suppose that at time zero the quantity a has a small but extensive spontaneous fluctuation a®. 
This will occur with exponentially small probability, but when it occurs the subsequent evolution of a(t) will be 
independent of the thermal noise, i.e. a(t) = E(a(t)\ao), where we denoted by E(-\ao) the conditional expectation 
over the trajectories for fixed ao at time zero. As a$ is small, we can write 

a(t) = E(a(t)\a ) = A(t)a (20) 

Denoting by E ao (-) the average over the distribution of ao, and C' a (t) = E(a(t)a(0)) the correlation function, it 
follows that A{t) = C a (t)/E ao (0%). Notice that the typical values of ao entering in the correlation function are of the 
order y/N, while in the relation (|2^) we consider values of order N. The validity of the analysis above relies on the 
smoothness of the probability distribution of ao in the crossover region, assumption which is at the heart of linear 
response theory. 

As one is conditioning (^0|) by the value of ao only, then the overwhelming majority of the configurations C giving 
rise to the fluctuation are the ones "typical" of the restricted canonical distribution 

e-^8(a(C)~a ) 



JdCe-P H ( c )5(a(C)-a ) 
which is equivalent to 

„-/3(ir(C)- 7 a(C)) 



(22) 



in which 7 is fixed by: (a) 7 = ao- Since to the linear order in 7 we have (a) 7 = f3j(a 2 ) 7= o, it follows that the 
relaxation of a(t) induced by the field is given by: 

a(t) = faC a (t). (23) 

which is the fluctuation-dissipation theorem in its integral form. 

Here we would like to show how a generalized form of the regression principle holds in aging dynamics where the 
time scale separation suggests that, besides the fluctuations of the instantaneous magnetization M(u) one should also 
consider possible fluctuations of the running global magnetization m(u), defined as 



1 N 

m (*) = jV (24) 



i=l 



We consider the conditional expectation value of the magnetization at time t given small values of the instantaneous 
and running magnetizations M{u) and m{u): E (M (t)\M (it) , m(u)) . This can again be expanded to the first order: 

E{M(t)\M(u),m(u)) = A(t,u)[M(u) - m(u)] + B(t,u)m(u). (25) 

and the functions A and B can be fixed by a hypothesis of continuity, leading to 



E(M(t)\M(u),m(u)) 



C st{t~u) C ag {u/t) ' 
— [M(u) — m(u)\ H m(u) 

1 - qEA qEA 



(26) 



where we have used ((M(u) — m{u)) 2 ) = 1 and (m(u) 2 ) = SL ff^. 
Onsager's argument demonstrates two things: 

1. that the decay of a spontaneous fluctuation with time is governed by the correlation function. 

2. that a fluctuation induced by a conjugate field will decay as a spontaneous fluctuation if the probability distribu- 
tion defining the state of the system immediately after the induced field is turned off is equal to the unperturbed 
probability distribution projected on the hypersurface defined by the equations: 

N -, N 



i^m i («) = m(«); ^ £ - M(u), (27) 
1=1 1=1 

where now, m(u) represents the value of an average as (^|) for times immediately after the field is turned off. 
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This second condition is consistent with our scenario of quasi-equilibrium in the dynamical relaxation process. In the 
next sections we will deal with the problem of computing the slow part of the magnetization induced by a field. We 
will first discuss the case of SK-like models whose OTOs during the dynamical relaxation tend to the ground state 
values. Then we will discuss those systems where the asymptotic state is different from the ground state. In this case, 
the argument is further complicated by the extensive multiplicity of the threshold states. 



V. THE CASE OF SK-LIKE MODELS. 

We first recall that the equilibrium analysis of these models |l7j in the "one step replica symmetry breaking" 
approximation determines the multiplicity of states at low free-energy F fl§f| as 

Af(F)dF = e ^{F-Fas) dF (2 8 ) 

where Fqs is the ground state free-energy and x is the Parisi parameter in this approximation. 

We then quote from the dynamical solution the expression of the magnetization in the TRM experiment described 
in the previous section: 

M(t) = C st {t - u)/3h + C ag {u/t)l3hx. (29) 

The comparison of this with equation (|2^) tells us the following remarkable fact: the action of an external field h 
from time t w to u produces at u a state of the system (in the sense of a measure in the microscopic variables) which is 
identical to the one we can obtain through infinite realizations of the thermal noise and selection of those trajectories 
with 

M(u) - m(u) = /3h(l - q BA ) 

m(u) ~ fihxqEA- (30) 

Therefore, thanks to the use of Onsager's postulate it is enough to calculate the response at time u immediately 
after the magnetic field has been turned off. We have assumed t w and u sufficiently large so that the system is in a 
quasi-state with free-energy F slightly larger then the one of the ground state. Eqs. ([30|) separate the response to the 
magnetic field in two components: a) inside the same quasi-state the more probable configurations will change and 
b) the quasi-state will change. The response a) is the equilibrium intrastate response and is trivial. To isolate b) we 
imagine turning off the magnetic field at time u and then waiting a finite time At such that C s t(At) is qEA while still 
At/u is zero. Then we know that the system has gone from one quasi-state at time t w to another at time u + At both 
defined with zero magnetic field. The distribution of (zero magnetic field) quasi-states with this free-energy is given by 
(p8|). Each of them may have a magnetization, uncorrelated from the free-energy and with variance (to 2 ) = qea/N. 
The typical number of quasi-states with free-energy density F and magnetization m is therefore given by 

Af(F, to) = e 0x{F - FGs) e- N ^ (31) 

implying that 

m 2 N 

S(F,m)=px(F-F GS )-- >0 (32) 

iqEA 

We first note that if we send u to infinity before sending h to zero, i.e. we consider fields such that the induced 
magnetization to verifies f3x(F u — Fqs) << m 2 N/2qEA we can derive the result ( |30|) in a quite straightforward way. 
In fact, we obtain that a non zero magnetization has to be compensated by an increase of free-energy so as to keep 
the configurational entropy S (F, to) non negative 

F-Fas = ^L. (33) 
2q EA px 

Along this line in the F, to plane the state with lowest total free-energy F — hmN has 

F = F GS + ?^^. (34) 

implying to = /3xhqEA- 
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The interpretation of this result is particularly illuminating. Turning on the magnetic field is a way of making 
energy available to the system. The thermal bath would normally absorb part of this energy. However this is possible 
only if the entropy of the system decreases in the process. This cannot happen here as by hypothesis the available 
entropy is much smaller than the one required to increase m. We conclude that the equilibration must occur only 
between the magnetic free-energy hmN and the unperturbed, zero magnetic field F. 

With this argument in mind we can now understand the limit more relevant to the dynamical approach. In this 
case F(u) — Fgs is extensive and large with respect to the potential energy introduced by the external magnetic field. 
In this situation there is, formally, enough entropy to allow the magnetization to reach the value of equilibrium with 
the thermal bath m = (3hqEA- However, with the same token one would argue that the thermal bath could have 
absorbed that entropy to decrease the spin-spin interaction energy. We know that this is not the case, or rather that 
entropy/heat is absorbed at a certain rate basically determined by the barriers. The external force is uncorrelated 
with the direction of relaxation of the system, and therefore it is reasonable to assume that the turning on of the 
magnetic field will not modify the rate of entropy decrease (heat transfer to the thermal bath). We conclude as before 
that the equilibration must occur between the magnetic free-energy hmN and F . In formulas if we call F h (u), F(u) 
the free-energy (associated with the inter spins couplings) that the system would reach in the presence of the magnetic 
field or in its absence at time u, then: 



so that: 



S(F h (u), m(u)) = S(F(u)) (35) 

(3x(F h (u) - F GS )-^— = Px{F{u) - F GS ) (36) 
The previous argument now follows minimizing F h {u) — Nhm(u). 

We remark that both entropy reductions refer to the same degrees of freedom, and therefore respond on the same 
time scale. The result is that the thermal bath acts as if it was uncoupled while the two forms of (free-) energy 
mutually equilibrate.^] In other words the transition time to higher free-energy states is much smaller than the one 
required to go to equal or lower free-energy states. 

This argument is so crucial to our picture that we feel it necessary to try to confirm it with a detailed model of the 
dynamical process. 

Let us imagine the dynamical trajectory from a (large) time u to a time t such that C(t, u) 0. We discretize the 
dynamics in k steps such that u = to < t\ < ... < tk = t such that C(tj+i,ii) = qEA — e with e small. At time U 
the system will have free-energy Fi and -F,+i — Fi will be small but extensive. The model we make of the dynamical 
process consists in assuming that when going from time ij to time U + \ the system can access different quasi-states 
with the lower one at free-energy i^+x and the higher ones distributed exponentially 

Af(F) = Q^( F ~ Fi +^ (37) 

while the probability of transition to a state with free-energy F is proportional to 

e-? F (38) 

The model is consistent for A < /3, (otherwise the free-energy would grow with time) and incorporates the following 
two features: 

1. The decrease in extensive free-energy is deterministic 

2. If we fix the initial condition, the increase in entropy in a single step is finite.^ In fact this can be calculated 
following the lines in for the Random Energy Model, with the result 



3 This represents an instance of the recent proposal that a system and a thermometer responding in the same time scale 
will equalize their effective temperatures fid]. In fact the inverse temperature of the magnetic field interaction energy is 
dS/dEh = d(—m 2 /2qEA)/d(—mh) — m/(qAEh) — fix. However our entropic interpretation suggests that time scales will 
strongly depend on f3x, the lower the effective temperature, the slower the evolution of the system; if two different aging 
systems starting with different effective temperatures and equal time scales are put in contact, they will quickly develop 
different time scales before equilibrating. 

4 The entropy about which we are talking here correspond to a the dynamical probability in which the initial condition is 
fixed, and is therefore increasing with time. 
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(39) 



In k steps the entropy generated will be kAS and therefore negligible with respect to Nm 2 / (2qEA)- Although non- 
extensive, kAS can be arbitrarily large, thus explaining the divergence of two cloned trajectories. In this model we 
have heavily used the self-averaging character of the macroscopic quantities along the trajectories. It is again clear 
that the only way to develop a magnetization is by compensating it with an increase in the (zero magnetic field) 
free-energy. 



VI. P-SPIN LIKE MODELS 



For p-spin like models even the limit t w — ► oo before h — ► is non-trivial. In fact it is well known that for this 
kind of systems the properties of the quasi-states encountered in the dynamics are closer and closer to these of the 
threshold TAP states, which have extensive configurational entropy Sth- If this entropy would be accessible in the 
dynamical process the equality (|l^) would be valid with = Sth- The condition that the total configurational 
entropy at time u be positive would then read 

m 2 N 

S th + (3x(F u -F th )-- >0 (40) 

which could be satisfied even if F u — F t h is small and negligible in front of — f^j ■ In more physical terms we can say 

that among the e th states there are e th ~ 2 iEA states with magnetization m. If all these states were available to a 
single trajectory the response would be normal, m = flhqEA- If we want the response to be anomalous we must show 
that the system, while wandering in phase space has no access to the configurational entropy.^ 

The logarithm of the number of states in the vicinity of any given TAP state has been computed by Cavagna, 
Giardina and Parisi in ref. p3] in the case of the p-spin model. Below threshold all states are isolated; there are no 
states closer than ci given distance Qea — Qmax- Qmax 

is a level-dependent overlap which tends to qea at threshold. 
Right at the threshold, the logarithm of the number of states as a function of the distance (qEA — q) is 

NZ(q)<xN(q EA -q) 5 . (41) 

Let us again imagine a discretization of the dynamics in which at each step the system can jump a distance 
8 = (q — qEA)- Then after n steps the log of the number of accessible states would be at most of the order of n 5 5 .[] 
On the other hand the distance traveled will be: A = n 8 if all the steps are in the same direction and A = y/n 8 if 
the steps are uncorrelated. In both cases it is easy to see that if we take the limit 8 — > and n — ► 00 fixing A we find 
that log(Af)/N goes to zero (n S 5 — > 0). 

Notice that the argument is based on the scarcity of states in the vicinity of a given state. This should be a generic 
feature for p-spin like systems other then the p-spin model. 

Having eliminated the configurational entropy from the balance, the argument proceeds as in the case of SK-like 
models. 

Let us conclude by pointing out that threshold states with large magnetization (of order /3/m?ba) do exist, but are 
non-critical in the presence of the field. Therefore with probability one such states would be isolated and unreachable. 



VII. SUMMARY AND CONCLUSIONS 



The main point of our analysis has been to give an explanation of the anomalous response function. We have 
found the physical origin of the equality between the FDR and the growth rate of the configurational entropy close 



5 A moment of reflection reveals that otherwise the system, wandering in such a large space, would pass to lower lying states 
and relax below F t h- 

6 This estimate could correspond to a severe double counting, as one can realize applying the estimate to finite dimensional 
Brownian motion. In infinite dimensional problems we expect it to give essentially the correct result. But, in any case, we only 
need it as an upper bound in our argument. 







to the asymptotic state. The value of the anomalous response can be traced to the lack of available entropy when 
the system is close to the low lying states. Our interpretation clarifies the relation among equilibrium properties and 
off-equilibrium dynamics. For p-spin-like systems we have argued that the extensive configurational entropy of the 
threshold states does not play any thermodynamical role. We have seen that the classical Onsager's argument on 
the equivalence between the regression of a spontaneous, noise-caused, fluctuation of the magnetization and the one 
induced by an external field can be generalized to aging systems. 

Our analysis can be summarized by saying that in aging systems the rate of entropy decrease is a function of age 
and does not change due to small forces. Thus the balance is always between the value of the unperturbed free-energy 
and the one of the perturbation, without taking into account the thermal bath. We expect this conclusion to hold 
also in short range systems with aging. 

In spin glass materials, one time observables (OTO) equilibrate, and the picture we have developed relates to the 
structure of configuration space close to the ground state. In aging experiment of structural glasses on the other 
hand, OTOs are far from their asymptotic values. Still, one can observe quasi-scaling aging dynamics on two time 
observables. The structure of the phase space visited on this time scale can not be related to "true" asymptotic 
properties of the system. We would like to speculate that even here the fluctuation-dissipation ratio, which could be 
a slowly varying function of time, is related to the derivative of available phase space with free energy also varying 
along the dynamical path. This could be true even if the system would eventually reach equilibrium on a different 
time scale where FDT is asymptotically obeyed. 

Finally the case of multiple time sectors or multiple replica symmetry breakings will need trivial modifications. 
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VIII. APPENDIX 



The aim of this appendix is two fold. We first show that in the p-spin model the derivative of the configurational 
entropy of the saddles is continuous at the threshold. Then, we prove that above Td the paramagnetic state can be 
seen, to the first order in T — Td as a collection of quasi-states identifiable with the points of least gradient of the 
TAP free-energy. 

Let us start from the expression of the TAP free-energy for the p-spin model |21[| 



F T Ap[mi = y/qSi] = E Q q~ 



(3 (1 - p (1 - q) q- 1+ P - qP) log(l - q) 
4 2/3 



(42) 



where Eq is the angular part of the energy as a function of the angular variables 5,;. It is well known that while 
one can find stationary points of the angular part for all the values of E in the range | £7 1 > — -^gs- Conversely, at 
finite temperature one finds solutions for the radial parts only in the range —E t h > \E a \ > —Eqs- The overwhelming 
majority of these solution are free-energy minima. 

The stationary points of the angular part for —E t h < \Eq\ turn out to be saddles, with a number of unstable 
directions which depends on Eq. The number of stationary points as a function of Eq is given by |22| 



V 



p2 z 2 



(-1+p) Z 2 



1 (P Z > 

log(— , 



(43) 



where z is an auxiliary variable given by 



■1+p 



V Eq 2 — Eth'' 

-1+p 



(44) 



For the saddles E > E th = - (^/2 

which appears in the computation 
determinant. As suggested in pl| 
parabolic shape 




the formula become complex. This is due to the fact that the Hessian 

2[ has negative eigenvalues and one has to compute the absolute value of its 
this can be done just taking the real part of expression (fl3h, which gives the 
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£(|£ | < -E th ) = -El ^—L + - lo g (p - 1) (45) 

An explicit computation using this formula shows that the i5o-derivative of (J|^) and j45| ) is continuous at the threshold 
energy. 

Let us now pass to our second task. We would like to identify the quasistates close to threshold as points of minima 
of TAP gradient. Unfortunately we were not able to prove this directly in the aging regime at low temperature, for 
we do not know how to compute the dynamical entropy. We start then from the observation that for temperatures 
higher, but close to Td one observes slowing down of the dynamics with time scale separation which is less and less 
ambiguous as T — > Td- So we can define dynamic quasistates even above T c , where the role of a small but finite T — Td 
is similar to the role of a finite t w in the low temperature dynamics. We put these quasistates in relation with the 
TAP free-energy, supposing that they coincide with the points of least TAP gradient for fixed internal energy equal 
to the paramagnetic value —(3/2. These are saddle points of the angular part, while the radial part is an inflection 
point, i.e. we fix q in the value of the minimum of cLFtap / 'dq. By explicit computation from ( fl2| ) and ( pl5| ) we find 
that the total free-energy Ftap — TH(Eq) is equal to the paramagnetic free-energy —0/4 up to terms which are of 
the second order in T - T d . For instance for p = 3 one finds that F TAP - TE(E ) = -0/4 - 8y/2/3(T - T d ) 2 . 
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